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We present deep and wide V, I CCD photometry of Ursa Major I (UMa I) dwarf spheroidal galaxy (dSph) in Local Group. The images 
of the galaxy were taken by Subaru/Suprime-Cam wide field camera, covering a field of 34' x 27' located at the centre of the galaxy. 
Colour-magnitude diagram (CMD) of the UMa I dSph shows a steep and narrow red giant branch (RGB), blue and red horizontal 
branch (HB), and main sequence (MS) stars. A well-defined main sequence turn-off (MSTO) is found to be located at Yqmsto ~23.5 
mag. The distance modulus is derived as (m - M)o = 19.93 +0.1 (corresponding to a distance D= 96.8 ± 4 kpc) from the V-band 
magnitude of the horizontal branch (Vo,hb = 20.45 ± 0.02). The mean metallicity of the RGB stars is estimated by the V-I colour to 

5_( . [Fe/H] 2.0. The turn-off age estimated by overlaying the theoretical isochrones reveals that most of stars in the UMa I dSph are 

formed at very early epoch (~ 12Gyrs ago). The isopleth map of stellar number density of the UMa I dSph, based upon the resolved 
star counts of MS, RGB, HB stars as well as blue stragglers (BS), shows that the morphology of the UMa I dSph is quite irregular 
and distorted, suggesting that the galaxy is in a process of disruption. The very old and metal-poor nature of the stellar population 
implies that the star formation history of this newly discoverd faint dSph may have been different from other well-known 'classical' 
dSphs, which show significant stellar population of intermediate age. The stellar population of the UMa I dSph closely resembles 
that of Galactic old metal-poor globular cluster, but its size is typical of Galactic dSphs (r c =188 [pc], ri/2=300 [pc]), and the shape 
, of its spatial density contours suggests that it is undergoing tidal disruption. These characteristics of stellar population and spatial 

distribution of the faint galaxies help us to know how they formed and evolved, and give a hint to the nature of the building blocks of 
0^ ' hierarchical galaxy formation. 
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■ Milky Way. Modern cosmological models based on the Cold Belokurov et al. |2007l Irwin et al. |2007l Walsh et al. 2007). 

Dark Matter paradigm demonstrate the importance of hierarchi- The Ursa Major I (UMa I) dwarf spheroidal galaxy locating 

cal structure formation in all scales. Galaxies like the Milky Way a b ut 100 kpc from the Milky Way is the first dSph discovered 

and M 3 1 form as a part of a local overdensity in the primordial by SDSS (Willman et al. |2003] |. The absolute magnitude of the 

matter distribution via the agglomeration of numerous smaller UMa I dSph is M v ~ -6.75, which is about 8 times less lu- 

building blocks which independently could develop into dwarf m i n0 us than the faintest "classical" Galactic dSphs, such as the 

galaxies if they could have been escaped from any take over by Draco, Ursa Minor, and Sextan s dSph . Measuring veloci ties of 

bigger ones. The relatively gas-rich dwarf irregulars still exhibit the UMa I stars, Kleyna et al. d2003T > and Martin et al. d2007b 

on-going star formation, while the dwarf spheroidals, being de- suggested that the UMa I dSph might be the most dark matter 

void of significant amounts of gas and dust, are now quiescent dominated object in the Universe. 

and are therefore, in principle, much simpler systems to study. xhe propeities of new ly discovered faint dSphs are poorly 

The proximity of Galactic dSphs offers a unique opportunity for under stood. Because of their faint luminosity and the size in the 

investigating galaxy formation and evolution in unprecedented sky ^ the stellar ages ^ me tallicities and structures of such dSphs 

detail by studying the photometric and spectroscopic properties are stm unclear from previous liming observations. The super- 

of stellar populations. wide field of yiew and en0U g n depth of the Subaru/Suprime-Cam 

offers a unique opportunity for successfully and efficiently ad- 



1. Introduction Analyses of resolvable stellar distribution in the Sloan 

Digital Sky Survey (SDSS) data archive have led to recent dis- 
The origin of Galactic dwarf spheroidal (dSph) galaxies should co veries of numerous very faint G alactic satelli tes (Willman et 
closely be related to the formation and evolution history of the a i. [2005] Belokurov et al. 120061 Zucker et al. I2006al |2006b; 



Send offprint requests to: S. Okamoto, dressing the photometric aspect of this research, 

e-mail: sakurako . okamoto@nao . ac . jp This paper is organised as follows: In Section 2, we describe 

* Based on data collected at Subaru Telescope, which is operated by the observation and data analysis procedure. In Section 3, we 

the Nationa Astronomical Observatory of Japan. present colour-magnitude diagrams (CMDs) of the UMa I dSph, 
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Fig. 2. Completeness of I (solid line) and V (doted line) photom- 
etry. The horizontal doted line represents the 90% level. 



Fig. 1. The V-band Suprime-Cam image of Ursa Major I dSph 
covering a field of 34'x 27'. North is to the top and east is to the 
left. 

and discuss the stellar population of the UMa I dSph. In Section 
4, we estimate the age and metallicity of the dominant stellar 
population. In Section 5, we analyse projected density profile 
based upon the resolved star counts. Then, we discuss the struc- 
ture and evolution of the UMa I dSph and summarise this study 
in Section 6. 



2. Observation and Data Reduction 

We obtained deep V, I CCD images of the central part of the 
UMa I dSph with the Suprime-Cam (Miyazaki et al. l2002b on the 
Subaru Telescope during nights of 2005 December 30 to 2006 
January 2 (PI. N. Arimoto). The Suprime-Cam consists of a 5 
x 2 arrays of 2048 x 4096 CCD detectors and provides a field- 
of-view of 34'x 27'wifh a pixel scale of 0.202". The nights of 
the observing run were photometric and the seeing ranged from 
0.5"to 0.8". To avoid saturation of bright stars, we took long and 
short exposure images with Johnson V (10 x 50s and 6 x 10s) 
and Cousins I (10 x 1 10s and 3 x 30s) filters. 

Data were processed using a pipeline software SDFRED 
dedicated to the Suprime-Cam (Yagi et al. 2002; Ouchi et 
al. [2004). Each image was bias-subtracted and trimmed, flat- 
fielded, distortion and atmospheric dispersion corrected, sky- 
subtracted, and combined in usual manner. Astrometric calibra- 
tion of each passband was based on a general zenithal polyno- 
mial projection derived from astrometoric standard star selected 
from online USNO catalog Figure Q] shows the V-band com- 
bined image of the UMa I dSph. Because of its faint luminosity, 
we cannot see the shape of the galaxy directly. 

For these processed images, the DAOPHOT in IRAF pack- 
age was used to obtain the point-spread-function (PSF) photom- 
etry of the resolved stars (Stetson |19871 l. In order to separate 
the point sources from the extended ones and noise-like objects, 
we used DAOPHOT parameter^ 2 and SHARP with the follow- 
ing cuts: x 2 < 4 and -0.4 < SHARP < 0.4. The position of de- 
tected stellar objects in each band was cross-correlated (within 
1") in order to make catalogues of long and short exposure im- 

1 http ://f tp . nof s . navy, mil/data/f chpix/ 



ages. Then, combining bright stellar objects (V < 21) detected in 
short exposure image and faint stellar objects (V > 21) detected 
in long exposure image, we obtained a candidate list of about 
2800 stars. 

Instrumental magnitudes of sources in the CCD images were 
calibrated to the standard Johnson-Cousins photometric system 
using the photometric standard stars of Landolt ( 1992), observed 
during the same nights. To confirm our photometric calibra- 
tion, we checked the consistency of the magnitude of stars com- 
monly detected in both long and short exposure images (typi- 
cally 20-22 mag stars). The difference between the magnitude 
of the long and short exposure is h.M S hort-u>ng < 0.02. The av- 
erage extinction, E(B-V)= 0.019 + 0.005 (corresponding to 
A y = 0.06 + 0.01, A 7 = 0.04 + O.Ol B in th e direction of the 
UMa I dSph is taken by Schlegel et al. ([T9981 . 

We derived completeness and photometric errors of our V, I 
photometry using the artificial-star test with the ADDSTAR rou- 
tine in DAOPHOT. We added 7000 artificial stars to each CCD 
image in each 0.5 magnitude interval from 18 mag to 24 mag 
and in every 0.2 magnitude interval from 24 mag to 26 mag. We 
processed the resulting images containing artificial stars in the 
same way as for the original images. The detection ratio of the 
test, N(recovered)/N(added), is plotted in figure [2] which shows 
that our photmetry is 90% complete at V=24.9 mag and 1=24.1 
mag, respectively, independently on position in the galaxy. 

The mean errors in the photometry are based on the differ- 
ence between the input magnitude and the output magnitude of 
the artificial stars. These errors are ploted in figure [3]in the next 
section. 



3. Stellar Population of the UMa I dSph 

Figure [3] show the de-reddened CMDs of the star-like objects 
found in the UMa I field. Error bars in both figures show the 
photometric errors at each magnitude level based on the artifi- 
cial star test. These data are the first to reach ~1 mag below the 
main sequence turn-off (MSTO), and one can easily see main 
sequence (MS) stars of the UMa I dSph. The major features of 
CMDs are as follows; (i) The red giant branch (RGB) is steep 
and narrow. The tip of RGB is buried in the foreground stars. 

2 The assumed extinction law is the R, =3. 1 (Cardelli et al. 1989 1 and 
the relation of A,/A v =0.594 (Schlegel et al. [E?98l 
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Fig. 3. Colour-magnitude diagrams of the UMa I dSph. Error bars show the photometric error at each magnitude level based on the 
artificial star test. The boundaries marked in the top-left panel are used to select stars belonging to the main sequence, giant branch 
and horizontal branch and blue strugglers of the UMa I dSph. 



The intrinsic width of RGB indicate a homogeneous chemical 
composition of the UMa I dSph. (ii) The horizontal branch (HB) 
stars are seen at Vo ~ 20.45 mag. (iii) Well-defined MS stars 
can be traced below Vo ~ 23 mag, with the MSTO located at 
Vo ~ 23.5 mag and (V— I)o ~ 0.5 mag. (iv) There are a few blue 
straggler (BS) candidates. Although it has not been addressed 
whether these stars are young MS stars or BSs, the paucity of 
BS candidates in UMa I dSph suggests that the galaxy had no 
multiple episodes of star formation. 

The CMDs are similar to those of metal-poor Galactic glob- 
ular clusters, implying that the UMa I dSph has very old and 
metal-poor stellar population (see Section 4). 

These CMDs include foreground stars and background ob- 
jects. The main comtamination is the Galactic foreground stars 
visible at (V-I)o > 0.6 mag. The contamination of the RGB re- 
gion in the CMD could in principle be high. However, the rela- 
tively high galactic latidude of the UMa I dSph helps in reducing 
this effect and the MS, HB features that we consider in figure[3] 
consist predominantly of stars belonging to the UMa I dSph. In 
our radial profile fitting in Section 5, we use the TRILEGAL0 
Galaxy model code (Girardi et al. 120051 ) to estimate the fore- 
ground star contamination. 

4. Distance, Age and Metallicity 

The RGB of the UMa I dSph in figure [3] is very thin and the 
tip of RGB can hardly be identified, so that the RGB tip of the 
UMa I dSph cannot be used as a distance indicator. We therefore 
have estimated the distance to the UMa I dSph using the V-band 
magnitude of the HB stars. We estimate the HB magnitude to 



3 http://trilegal.ster.kuleuven.be/cgi-bin/trilegal 



Vq.hb = 20.45+0.02. Since the UMa I dSph seems to have stellar 
population like Galactic metal-poor globular clusters, we assume 
theoretical value of absolute magnitude M^hb — 0.515*q5|, cor- 
responding to [Fe/H]= -2.0+U5 (Catelan et al. 120041) . which give 
the distance modulus (distance) as (m - M)q = 19.93 ±0.1 
(96.8 + 4 kpc). 

In the left panel of figure [4] we compare the V-(V-I) dia- 
gram with the fiducial sequences of metal-poor Galactic globular 
clusters M3, M5 and M92, taken from Johnson et al. ( 11998) : the 
metallicities of these clusters are [Fe/H] =-1.57, -1.99, -2.24, 
respectively (Harris 11996) . We adopt the distance moduli to 
(m-M)o = 14.99, 14.46, 14.67, and reddening corrections to 
E(B-V) = 0.0 1, 0.03, 0.02 for M3, M5 and M92, respectively 
(Harris [19961 The mean RGB and MSTO of the UMa I dSph 
are similar to those of M5 and M92, suggesting that the metal- 
licity of the UMa I dSph is [Fe/H] =* -2 and the age is as old as 
M5 and M92. 

Our metallicity estimate is consistent with previous photo- 
metric and spectroscopic studies (Willman et al. 120051 Martin 
et al. 12007) . and consistent with the metallicity we assumed for 
estimating the distance ([Fe/H]= -2.0*9.1). 

By overlaying theoretical isochrones on the deep CMD, we 
estimate the age of stellar population in the UMa I dSph. The 
isochrones in the right panel of figure [4] are interpolated Padova 
isochrones corresponding to metallicity [Fe/H] =-2.0 and age of 
10.0, 12.6, 14.1 Gyr (Girardi et al. 120021 . These isochrones are 
made by the interpolation with the Padova isochrones of metal- 
licity Z=0.0001 and Z=0.0004. The right panel of figure|4]shows 
that the fiducial sequence of MSTO of the UMa I dSph is best 
reproduced by the isochrone of 12.6 Gyr, consistent with the age 
roughly estimated by the comparison with Galactic old metal- 
poor globular clusters in the left panel. Consequently, most of 



4 



S. Okamoto et al.: The stellar population of Ursa Major I dSph 



> 




(V-I)o (V-I)q 



Fig. 4. Left: V-band colour-magnitude diagram of the UMa I dSph with the fiducial sequence of Galactic globular clusters (M3, M5, 
M92). Right: The theoretical isochrones overlaid in the V-band colour-magnitude diagram of the UMa I dSph. The lines stand for 
the interpolated Padova isochrones for [Fe/H]=-2.0 and age of 10.0, 12.6, 14.1 Gyr which were sifted according to the distance of 
the UMa I dSph. 



stars in the UMa I dSph are estimated to be formed at least 12 
Gyrs ago. 

The UMa I dSph has old metal-poor stellar population with 
the narrow dispersion of age and metallicity. This characteristic 
of the stellar population of the UMa I dSph seems to be different 
from "classical" dSphs, which show evidences of multiple stellar 
populations of various age and metallicity (e.g. Shetrone et al. 
I200T1 Ikuta & ArimotoEOOU Tolstoy et al. l2§04l . The UMa I 
dSph is well-described by a single stellar population, as a typical 
globular cluster. 

5. Structural properties of the UMa I dSph 

Even though its stellar population closely resembles that of 
Galactic old metal-poor globular clusters, the spatial distribu- 
tion of the UMa I dSph shows different properties in morphol- 
ogy and size. Figure|5]shows the spatial distribution and the den- 
sity contour map of stellar sources, which were selected as can- 
didate stars of MS, RGB, HB and BS from the CMD shown 
in figure [3] using the selection boxes marked. These stars are 
binned in 30"x30"boxes and smoothed by the Gaussian kernel 
of bandwidth=0.6'to draw the lower panel of figure[5] The plot- 
ted contour levels are 1-10 cr above the background level. The 
isodensity contours covering the central 730 pc x 670 pc of the 
UMa I dSph, have quite irregular shapes, stretched toward the 
upper left and lower right on figure [5] It corresponds to the tan- 
gential direction toward the Galactic centre. 

We also compare the spatial distributions of subgiant stars 
(23.0 < V < 23.4) and faint MS stars (24.0 < V < 24.4) to study 
the mass segregation in the UMa I dSph, but there is no differ- 
ence in both stars. 



To calculate the structural parameters, we use the candidate 
members in figure [3] to derive the centroid from the density- 
weighted first moment of their spatial distribution, and average 
ellipticity and position angle using the three density-weighted 
second moments (e.g. Stobie I1980I I. The radial profile shown 
in figure [6] is derived by calculating the average number den- 
sity within elliptical annuli. The number density of foreground 
stars in the direction of the UMa I dSph is estimated as 0.05 
arcmirT 2 by using TRILEGAL (Girardi et al. 120051 We fit 
the radial profile with standard exponential and Plummer mod- 
els, and estimate the half-light radius (ri/2) of both models as 
11.30'+ 0.50'and 10.42'+ 0.76', respectively. Hereafter the ra- 
dius is the elliptical radius r — [x 2 + y 2 /(l - e) 2 ] 1 ^ 2 , where e is 
the ellipticity of the galaxy, x and y is the coordinate aligned with 
the major and minor axis, respectively. The estimated half-light 
radius is larger than previous estimation of ri/2 ~7.75'(Willman 
et al. 120051 1. probably thanks to our much deeper images. The 
field level of Willman I (Willman et al. 2006) which lies within 
few degree from the UMa I is within the error bars of field level 
in our profile. 

The half-light radius of the UMa I dSph corresponds to 300 
pc, which is an order of magnitude larger than that of the largest 
Galactic globular clusters. The total luminosity of the candidate 
members is used to estimate the lower limit of the central surface 
brightness to fioy = 29.5 + 0.2 [mag/arcsec 2 ]. These best-fitting 
structural parameters are listed in table Q] 

The structural parameters allow us to constrain the dark 
matter content of the UMa I dSph, assuming that it is in 
virial equilibrium. The mass-to-light ratio of a simple stel- 
lar system with symmetric distribution can be estimated as 
M/L=977<x 2 llnGrhbLo, where 77 is a dimensionless parameter de- 
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Fig. 6. Raidal profile derived by calculating the average num- 
ber density within elliptical annuli. The effect of foreground 
stars is corrected by using Galactic model (Girardi et al. 
2005 ). Best-fitting exponential profile with the effective radius 
r e =6.73'+0.30'and Plummer profile with the Plummer radius 
b=10.42'±0.76'have been overlaid as dashed line and doted line, 
respectively. 

Table 1. Structural parameters of Ursa Major I dSph 



Parameter 


Value 


Coordinates (J2000) 


10h34m44.4s, +51°55'33.9" 


Position angle 


78°. 12 


Ellipticity 


0.54 


ri/2 (exponential) 


11.30 ±0.50 [arcmin] 


1 1/2 (Plummer) 


10.42 ± 0.76 [arcmin] 


floy (exponential) 


29.5 ± 0.2 [mag/arcsec 2 ] 


yUo.y (Plummer) 


29.7 ± 0.4 [mag/arcsec 2 ] 


(m-M)o 


19.93 ±0.1 


Distance 


96.8 ± 4 [kpc] 



Fig. 5. Upper: The spatial distribution of the stars selected by 
the dashed lines in figure [3] which coverd the central 730 pc x 
670 pc of the UMa I dSph. Lower: The isodensity contour of 
the CMD-selected sources. The plotted contour levels are 1-10 
<x above the background level. 

pendent on the luminosity distribution, <x is a central velocity 
dispersion, r/ t b is the half-light radius, and So is the central sur- 
face brightness (Richstone & Tremaine |19861 l. Using the veloc- 
ity dispersion values <x = 9.3[km/s] (Kleyna et al. 2005) and 
cr = ll.9[km/s] (Martin et al. 120071 1, the mass-to-V band light 
ratio of the UMa I dSph is M/L ~ 1500 and 2500M Q /Z<r> re- 
spectively. This estimate implies that the UMa I dSph is a highly 
dark matter dominated galaxy. However, the shape of the spa- 
tial density contours suggests that the UMa I dSph is undergoing 
tidal disruption. If this is the case, the galaxy may not be in virial 
equilibrium. 



6. Discussion and summary 

We have used Subaru/Suprime-Cam to obtain multi-colour 
imaging of the UMa I dSph which is sensitive enough to de- 
rive the stellar population below the MSTO and wide enough 
to study the spatial distribution of stars in the galaxy. Our deep 
and wide photometry allows us to study more accurately than 



previous study. The distance modulus of the UMa I is derived as 
(m-M)o = 19.93+0.1, corresponding to 96.8+4 kpc. We identify 
a number of stars in several different evolutionally phases, MS, 
RGB, HB and BS, which are characterized as a metal-poor sin- 
gle old population. The mean metallicity is estimated by the V-I 
color of RGB stars to [Fe/H] — 2.0, and the turn-off age reveals 
that the most of stars in the UMa I dSph are formed at least 12 
Gyrs ago. 

Most "classical" Galactic dSphs have a sign of multiple stel- 
lar populations - different spatial distributions of red and blue 
HB stars (e.g. Harbeck et al. 120011 Tolstoy et al. [2004b , the 
[ff/Fe] ratio of most stars so far observed in the dSphs lower 
than Galactic halo stars of similar metallicity (Shetrone et al. 
1200 U 120031 ), which implies that star formation in such dSphs 
lasted for several Gyrs (Ikuta & Arimoto 2002). On the other 
hand, the UMa I dSph seems to have a metal-poor single old 
population, and some of the other newly discovered dSphs also 
show a sign of similar stellar population (Bootes I: Belokurov 
et al. 120061 Coma Berenices: Belokurov et al. 120071 Hercules: 
Coleman et al. 120071 ) and appear to have a different property 
of stellar population from well-known "classical" dSphs. Why 
these dSphs are so faint and have single old population? If the 
SN-driven winds blow the gas earlier in such a low-mass galaxy 
because of the relatively shallow potential, the stellar population 
may keep the information of the initial star formation activity of 
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the galaxy. The abundance patterns of these galaxies allow us to 
reveal the nature of galaxy evolution. 

Although the CMD of the UMa I dSph shows a single epoch 
of star formation as a metal-poor globular cluster, its spatial size 
is typical of Galactic dSphs (ri/2=300 [pc]), and the shape sug- 
gests that it is undergoing tidal disruption. The morphology of 
UMa I dSph is quite irregular, stretched and distorted. Similar 
feature is found in the Hercules dSph which shows extremely 
flat structure (e ~ 0.65) (Coleman et al. 120071 1. If the tidal force 
could lead such an elongated shape, these galaxies have a high 
orbital ellipticities and underwent tidal disruption at the last peri- 
center passage to the Milky Way. 

With the usual assumptions of symmetry and virial equilib- 
rium, the UMa I dSph seems to have high M/L and would be 
the dark matter dominated object. If so, the UMa I dSph is pre- 
sumed to be a long-lived object with distorted morphology. It 
is unclear whether the UMa I dSph have been caught in the act 
of dissolving into the Galactic halo recently, or maintains such 
elongated shape with dominated dark matter for a long time. The 
high-dispersion spectroscopic data will necessary to illuminate 
the chemical and dynamical evolution of the galaxy. 
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